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Abstract 

Within the framework of a generic generally covariant quantum theory we introduce the loga- 
rithmic correction to the quantum wave equation. We demonstrate the emergence of the evolution 
time from the group of automorphisms of the von Neumann algebra governed by this non-linear cor- 
rection. It turns out that such time parametrization is essentially energy-dependent and becomes 
global only asymptotically - when the energies get very small comparing to the effective quantum 
gravity scale. We show how the logarithmic non-linearity deforms the vacuum wave dispersion 
relations and explains certain features of the astrophysical data coming from recent observations 
of high-energy cosmic rays. In general, the estimates imply that ceteris paribus the particles with 
higher energy propagate slower than those with lower one, therefore, for a high-energy particle the 
mean free path, lifetime in a high-energy state and, therefore, travel distance from the source can 
be significantly larger than one would expect from the conventional theory. Apart from this, we 
discuss also the possibility and conditions of the transluminal phenomena in vacuum such as the 
"luminal boom" and Cherenkov-type shock waves. 
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I. INTRODUCTION 



In the conventional quantum mechanics the linearity of the wave equation is something 
which is implicitly presupposed, yet the possibility of the non-linear generalization has not 
been ruled out by experiment |l| . From the theoretical point of view, there exist arguments 
that a nonlinearity in general can lead to the violation of locality via the Einstein- Podolsky- 
Rosen (EPR) apparatus Q]. However, it has been also pointed out that instead of this 
nonlinearities can lead to communications between branches of the wave function, they can 
)e large in a fundamental theory yet be unobservably small when measured experimentally 
3)]. Afterwards it was also shown that the locality is not violated for a large class of 
nonlinear generalizations, including the one which is being discussed here, see Ref. (4], 
and references therein. In any case, the linearity requirement becomes a rather strict and 
unnecessary assumption if one expects quantum mechanics to be valid on a wide range of 
scales: for instance, the modern theory of quantum gravity is believed to be essentially non- 
linear - because the propagating particle will cause the quantum fluctuations in gravitational 
medium which will react back. 

Some non-linear extensions of QM have been already proposed - for instance, in Ref. [^j 
authors studied a family of the non-linear wave equations for non-relativistic QM associated 
with unitary group of certain non-linear gauge transformations of third kind - those which 
leave the positional probability density invariant. Another general approach to including 
non-linearity into QM is based on generalizing the quantum phase space to the class of 
Kahler manifolds which admit certain Hamiltonian flow {(]]. 

On the other hand, those who want to add non-linear terms into the wave equation 
inevitably arrive at the problem of choice: it seems that there exists a huge variety of the non- 
linear corrections which can be added without undermining the pillars of the conventional 
QM - the concepts of the physical state, probability, observables and measurement Q. 
Thus, the necessity of non-linear corrections encounters the practicality issue: what are the 
before unsolved problems of the conventional QM which can be cleared up by introducing 
non-linearity? 

In present paper we do not ab initio postulate the linearity of quantum wave equation 
at all energy scales. Instead we consider the including of one particular nonlinear term in 
a generally covariant way. We choose the logarithmic one for in the flat-spacetime limit 
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We show that this 



it would not induce correlations for non-interacting systems 
nonlinearity naturally transforms into the evolution time derivative which becomes global 
and energy-independent only in the low-energy limit - when the energies are small comparing 
to the effective quantum gravity scale. Then we discuss the vast observational implications 
of our model. 

II. NON-LINEARITY AND ADDITIVITY 

The formal structure of a generally covariant quantum theory is as follows. Let T ex 
be the space of solutions of the generally covariant equations of motion endowed with a 
degenerate symplectic structure defined by these equations. The degenerate directions of 
this symplectic structure integrate in orbits and the solutions which belong to the same 
orbit must be physically identified. The orbits form a symplectic space T - a fully covariant 
object which becomes the physical phase space of the theory. The set A = C°°(r) of the 
real smooth functions on T, called the physical observables, form an Abelian multiplicative 
algebra. These observables are regarded as classical limits of the non-commuting quantum 
observables whose ensemble forms the non- Abelian C*-algebra A. Since A is a non- Abelian 
algebra under the Poisson bracket operation one can assume that A is a deformation of a 
subalgebra of the classical Poisson algebra. Note that in A in general there is no defined 
Hamiltonian evolution or representation of the Poincare group, therefore, the time evolution 
is only determined by the dependence of the observables on clock times. Suppose, at the 
classical level the evolution is governed by the constraint (more precisely, a combination 
of spacetime diffeomorphisms' constraints) % ~ which vanishes in a weak sense. In the 
quantum case one must define first a state u: A — > C, a positive, linear and normed 
functional on A. From the Gelfand-Naimark- Segal theorem it follows that there exists 
exactly one representation n^: A — > End% w of the algebra A on a Hilbert space "H^, and 
the vector £ w G W. w such that: (i) lAn(-K w (A)t, w ) = U^, (ii) u(a) = (^(a)^, ^) for every 
a E A [9|. Then the evolution of physical states is governed by an appropriately chosen 
operator %. 

To begin with, by G we denote the wave functional which describes the state of 
the dynamical system. Then in the generalized Schrodinger picture the quantum evolution 
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equation can be written in some representation as (we consider pure states for simplicity): 



V = 0, 



(1) 



where the first term in brackets is essentially the above-mentioned combination of constraints 
% quantized as in the conventional formalism. Its explicit form is determined by a concrete 
physical setup and thus will not be important for us here - we simply assume that it can be 
consistently defined. This will make our following results largely model-independent from 
spacetime-formulated theories. The other term, F(^f), is not present in the conventional 
quantization procedure. As to preserve the probabilistic interpretation of ^ (the physical 
states are actually not vectors but rays), we assume that F depends not on \1/ alone but 
rather on its complex square. Notice also that this term does not interfere with % as it 
describes the self-interaction of the wave functional and thus is inherent only to the way we 
define the quantum wave equation. Therefore, we write Eq. (CO as 



n + F{ P ) 



* = o, 



(2) 



where p = \^\ 2 - What is the explicit form of the operator Fl 

Suppose the system described by \1/ consists of two separated distinguishable subsystems, 
described by wave functionals and ^2, respectively, which obey the wave equations 



Ui + F( Pi ) 
form 



0. 



1,2, where we denoted pi = \^i\ 2 . We assume that F has the 



F(p) = -/3- i ln(fip) 



(3) 



where (3 and fl are arbitrary positive constants with the dimensionality of the inverse energy 
and space volume, respectively. This function is a general solution of the following algebraic 
equation 

F(p l p 2 )=F(p 1 )+F(p 2 ), (4) 

if has the same value for both subsystems, therefore, for two uncorrelated subsystems, 
when the wave functional of a whole system becomes the product 1^) = <E> ^2), the 
overall quantum wave equation turns to: "Hi 2 + -F(pi) + F(p 2 ) (^1^2) = 0, thus the non- 
linear part in this special case becomes a plain sum of the subsystems' ones. Therefore, the 
non-linear term obeying Eq. (T4J) introduces on its own no additional correlations between 
the uncorrelated subsystems for which (3 takes the same value. Later it will be shown that, 
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for instance, in the weak-gravity limit the relative difference between any two /3's is given 
by A/3//3 ~ AE/Eqq, where Eqq is certain very large energy scale and AE <C Eqg is 
the difference of the subsystems' energies, and therefore, in that limit Eq. (J3D holds with a 
high degree of precision. Thus, the logarithmic non-linearity on its own does not break the 
energy additivity and separability of non-interacting subsystems in non-relativistic quantum 
mechanics. Therefore, from the viewpoint of the correspondence principle the logarithmic 
term is safe to include into the quantum wave equation. 

Gathering all together, from Eqs. (j2J) and in the position representation we obtain 
the following quantum wave equation: 



which is the quantized version of the constraint H ~ in generally covariant theories. This 
equation can be formally written as 



space but the latter is not a problem because even in the conventional quantum formalism 
the space of physical states is already a projective Hilbert space (rays), due to the normal- 
ization constraint. Of course, to preserve interpretation of \I/ as a wave function, it must 
be normalizable and also the corresponding probability density must obey the conservation 
law which is the case for the logarithmic nonlinearity [8j. Otherwise, the nonlinear quantum 
wave equation can be interpreted only as an effective one, one example to be the Gross- 
Pitaevskii equation (yet the term "effective" should be treated with care here because in 
quantum gravity, unlike condensed matter, the background medium cannot be eliminated). 

To summarize, the nonlinear term ([3]) has a number of physically meaningful distinctive 
properties which make it a very probable candidate for a nonlinear correction to the quantum 
wave equation: (i) it would not break the general covariance of a theory because it does not 
interfere with the spacetime diffeomorphisms' constraints (yet, of course, it takes them into 




(5) 




(6) 



where is the Hermitian operator defined in the following way 10j: let us consider the 
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account - through the dependence on the wave functional which, in turn, is a solution of the 
wave equation), (ii) it does not violate locality via EPR-type correlations, as discussed above 
4|, (iii) if wave function satisfies Eq. (J5J) then this function times any constant is also a 
solution provided that %' is shifted by a constant, (iv) symmetry properties of wave functions 
with respect to permutations of the coordinates of identical particles are not affected by the 
nonlinear term, (v) in the non-relativistic limit the equations for the probability density and 
current are not altered, energy is additive and bounded from below, free-particle waves are 
localized. The non-linear corrections of other kind can also be included into Eq. but they 
would violate these properties, therefore, in what follows we focus only on the non-linear 
term Q. Below it will be shown that it has some additional properties which make it very 
useful in quantum gravity theories. 

III. MODULAR GROUP AND EVOLUTION TIME 

The physical meaning of the non-linear term in Eq. ([6]) becomes a bit more clear when 
we go to the non-relativistic limit. Its averaging in the position basis yields 

(/T 1 ^) = -/T 1 J \^\ 2 ln(tt\y\ 2 )d 3 x = T*S*, (7) 

where the functional is the Shannon-type entropy which can be formally assigned to a 
single quantum particle, = (ksP)~ l , is the Boltzmann constant. From the analysis of 
its discrete counterpart one can infer that reaches a minimum on the delta-like probability 
distribution (which corresponds to a classically localized particle) and maximizes on the 
uniform one [10|]. This is what the information entropy usually does, therefore, can be 
interpreted as a measure of the particle's "smearing" over space and corresponding quantum 
uncertainty. Such quantum-mechanical entropy is a purely information one and should not 
be confused with the von Neumann entropy - the latter would vanish for pure states and 
;hus neglect an inherent quantum- mechanical uncertainty of the outcome of a measurement 



ll| . Perhaps, the notion of entropy closest to S# would be the one proposed by Wehrl for 
the coherent states [121 ] . 

Further, strictly speaking, Eq. ([6]) contains no evolution time parameter. It is the 
well-known property of fully covariant theories that H contains the geometrical time-like 
coordinate at most but the dynamics cannot be formulated in terms of a single external 
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time parameter. This creates the conceptual difficulties because the full dynamics of the 
quantum gravity theory can not be consistently defined in terms of other available notions 



of time, such as the proper time, because they are state-dependent |13j . They can serve as 
the evolution times only for dynamical theories on a fixed spacetime geometry but not for 
the dynamical theory of geometry which the quantum gravity is supposed to be. 

One way to define the evolution time without invoking assumptions for spacetime geom- 
etry is to formulate it based on notions of a statistical and thermodynamical nature. To 
our knowledge, the initial idea was proposed in Ref. [l^] - to define the evolution time as 
the vector flow of the Gibbs state on the constraint surface. Is there any way to derive the 
notion of the evolution time from the quantum wave equation itself? We showed before that 
the nonlinear term from Eq. 03) can be interpreted as a kind of entropy, at least in the 
non-relativistic limit. Our next step will be to show that this term gives rise also to the 
evolution time in quantum theory. 

Luckily, the mathematical background for justifying this has been already developed, 
both for the conventional spacetime geometry and for the non- commutative one |is| . 
Recalling the notations above, let 72 be the von Neumann algebra generated by the rep- 
resentation of A, 7r w (^4), on the Hilbert space. Then the Tomita-Takesaki theorem asserts 
that the mappings of the von Neumann algebra on itself, a T : 72 — > 72 (r 6 M), of the 
form a T (b) = A~ tT 6A* r , b e 72, where A is a self-adjoint positive operator, form a one- 
parameter group of automorphisms of 7Z, called the modular group of the state u. This can 
be equivalently written as 

6=^-a T (6)| T=0 = i[6 J lnA], (8) 

(IT 

thus, the "modular" time a T can be already regarded as the time we are needed in but it 
seems still state-dependent so a bit of final tuning must be made. Some of automorphisms 
on 1Z are inner-equivalent, the set of their equivalence classes forms a group of outer auto- 
morphisms Out (7£) . Then the modular group a T of a state u projects down to a nontrivial 
group a T G Out (72). The state-independent characterization of time is then proven by the 
cocycle Radon-Nykodym theorem from which it follows that a T does not depend on choice 
of u in 72. 

Further, from Eq. fl6]) we obtain 

0=[b,H'] = [b,H]+p- 1 [b,SI 9 ], (9) 
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for any b E 1Z. Then, observing that the product e trn ' b e iTH ' is equivalent to e iTS *'Pbe ,rS */0 
on the constraints' surface, we can assume that 

A oc exp Si},, (10) 

from which, using Eq. (j8]), we obtain the equation of motion in the generalized Heisenberg 
picture: 

in) 

where the derivative is understood in a sense of the commutator from Eq. ([8]) with the 
generator given by Eq. ([TO]) . This essentially means that the dynamics described by the 
"stationary" operator M! is equivalent to the evolution governed by "H with respect to the 
evolution time j3r. In other words, the logarithmic nonlinearity can be "used up" for creating 
the time evolution of a generally covariant theory: theory containing such nonlinearity but 
without the observer- independent time evolution is equivalent to the linear theory with the 
evolution time defined by the modular group. As shown below, only in the low-energy limit 
this time parametrization becomes global. 

IV. DISPERSION RELATIONS AND OBSERVATIONAL TESTS 

From Eq. ([6]) one can see that the dispersion relation for a particle in vacuum is being 
deformed by the non-linear term. This is not surprising though as the quantum gravity is 



expected to give rise to such corrections 16Ml8| because the gravitational medium contains 



quantum fluctuations which respond differently to the propagation of particles of different 
energies - the phenomenon somewhere analogous to propagation through electromagnetic 



plasmas 19]. The full treatment of this problem is impossible without establishing a concrete 
quantum gravity model. Yet, some features are model- independent and can be clarified 
already at this stage. 

First, despite the functional form of the non-linear term is universal for all dynamical 
systems the constant is not a fundamental one hence depends on dynamical characteristics 
of a system. To find its physical meaning, we go to the flat-spacetime limit and consider a 
norm-preserving splitting of an arbitrary wave function into N non-overlapping parts of the 

N N 

same form as the initial wave function: — > ^/pl^{x — x{) where Ylpi = 1. Then 

i=i i=i 
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from the averaged Eq. ([6]) one obtains that the change of energy during the process amounts 



N 



to SE = —/3 1 Pi In pi, and thus is a measure of the binding energy, 

(3 oc 1/5E. (12) 

Second, from Eq. fill I) one can immediately see that for any two dynamical systems the 
relation 

{5\dr\ = (3 2 dr 2 , (13) 

must hold, thus, if two systems have different /3's then their evolution time scales must differ 
as well. 

Now, suppose that some two particles are products of the reactions happened inside a 
compact region of space. In the process these particles receive certain amounts of energy, 
E\ and E 2 , respectively. From the previous two equations we obtain that the ratio of their 
evolution time scales is given by 

at least in the leading-order approximation. Here Eq is the energy of vacuum of a theory, 
we imply that \Eq\ — Eqq where Eqq < 10 19 GeV is the effective quantum gravity energy 
scale. The large value of the latter explains why our (non-relativistic) notion of time is 
global and energy- independent: in the low-energy regime the characteristic energies of any 
two particles become small comparing to the quantum gravity energy scale, therefore, the 
difference between any two /3's also gets vanishingly small, ((3 2 — ~ (E 2 —Ei)/ Eqq, and 
when the value of /3 is essentially the same for all dynamical systems then this constant can 
be absorbed into the time parameter. Thus, the reason why the logarithmic nonlinearities are 
not observed in current quantum-mechanical experiments is not because of their smallness 
but because they act as the time derivatives with individual scale factors which are essentially 
indistinguishable in the low-energy regime. 

Another observation can be made if one recalls that the correspondence principle implies 
that for elementary particles their "modular" times must synchronize with their proper times 
in the weak-gravity limit. Then Eq. ( |T4l) immediately reveals the presence of the Lorentz 

invariance violation (LIV): it is not the conventional line element of spacetime which is 

□ 

invariant but the one multiplied by an energy- dependent function [20( . However, in the low- 
energy limit this circumstance does not bring any phenomenological difficulties because for a 



conformally flat spacetime (such as the Friedmann-Lemaitre- Robertson- Walker one) one can 
still define the representations of the Lorentz group in the vierbein basis, therefore, the main 



notions of particle physics are preserved [21|. In any case, LIV is an expected phenomenon 
in quantum gravity - the nontrivial vacuum creates a preferred frame of reference. 

Further, our particles fly off, travel across the space and eventually get caught by a remote 
detector. If they are of same kind, initially were emitted approximately simultaneously with 
equal velocities and their travel conditions were similar then what kind of differences between 
them are our detectors supposed to catch? Once again, the exact quantitative results are 
impossible without employing the concrete model of quantum gravity yet some heuristic 
analysis can be done in the flat-spacetime approximation (the cosmological corrections are 
not considered for now as to avoid certain confusions). 

Making in Eq. f fT4|) the transition from proper time to the distant observer coordinate 
time t (with cosmological effects that would be the comoving time), we obtain 



v x jc 2 -v\ E 2 - Eq 



(15) 



v 2 y c 2 - v\ E 1 - Eq 

where Vi = dx/dti. This equation reveals the following subtlety: if for our future purposes 
we assume that the particles are essentially relativistic or even ultrarelativistic and also that 
their velocities are nearly the same, then the value of a square root in the equation above 
crucially depends on whether the ratio V\jc tends to unity "stronger" than v\jv 2 - 

Thus, there exist two limit regimes of analysis: linear or standard relativistic - when 
the ratio Vi/v 2 approaches one "stronger" than V\jc does, and non-perturbative or extreme 
ultrarelativistic - when it is other way around. When analyzing concrete experimental data, 
one should look at these two ratios to decide which regime s/he is next to. Especially one 
should be careful when the cosmological effects are taken into account because these ratios 
may vary as the particles propagate. 



A. Linear regime 

In this case the square root in Eq. f|T5|) can be well approximated by one, therefore, under 
the above-mentioned assumptions we obtain: 

v 2 £i\ — J^o 
10 



thus, under equal conditions the particle with lower energy travels faster, so we can write 
v ~ v / (1—E/Eq) where = ve/e =q- The difference in their arrival times is proportional 
to their energy difference: 

For instance, for the photons we obtain At pa ^ AE, where L is the distance from a 
distant observer to the actual place where the reactions happened. This is what is being often 
observed about cosmic ray photons coming from the very distant Gamma-ray bursts (GRB), 
the highly energetic explosions of massive stars in galaxies: in the linear approximation the 
difference in the arrival times of photons is proportional to their energy difference. For 
instance, during the exceptionally luminous GRB 080916C, distant from us as far as 12.2 
billions light years, the first photons with energies above 1 GeV started to arrive only after 
ten seconds after the trigger, e.g., the 13.2 GeV photon had arrived after 16.54 s (the duration 



of the whole event itself was few tens seconds) [221 ] . The cosmological-scale remoteness of 
this and some other GRB's plays a crucial role here: from the last formula it is clear that 
L/c should be very large so as to win over the huge number Eqq/AE and produce an 
experimentally detectable effect. 

The observational predictions can also be formulated on language of the (deformed) 
dispersion relations for particles in vacuo. From Eq. ( fT6l) we obtain 

Av = v 2 - vi _ v 2 ^ 
AE E 2 — E\ Eq — E2 

therefore, if E <C \E \ ~ Eqq then dv/dE pa Av/AE ~ —^v/Eqq, where we assume 

£ = — sign (E ) = ±1. It means that the velocity can be written as an exponent of energy 

and thus is a linear function of E in the leading order: 

v/c ~ exp (-tE/Eqa) = 1 - + 0(E 2 /E 2 QG ), (19) 

-^QG 

and same result can be obtained if we look for a Taylor series solution of Eq. ( fl6l) : 

N 

v i/ c = S a n(Ei/E Q ) n , where N > 1 is an approximation order, a n are energy-independent 

n=0 

constants to be determined. The boundary conditions are determined by physics - for in- 
stance, for the case of photons they would be: Vi = c when Ei/Eo — > 0. 



With Eq. (fi~9l) in hands we can immediately recover the results of Refs. 18|, |23|, l24 |. 
Indeed, the resulting dispersion relation for the rotationally invariant case, BE /dp = v(E), 
integrates to 
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c 2 (p-Po) 2 -E 2 QG (l-^ E ^) 2 = E 2 



E 



1 + + <D(E 2 /E' QG ) 



QG 



(20) 



where po is an integration constant. In the limit po — > the relation reduces to the 
dispersion relation for massless particles in the effective quantum gravity theories based 
on K-deformations of Poincare symmetries in which the time coordinate does not com- 



mute with spatial ones 
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26]. On the other hand, with the appropriate choice of the 



integration constant the last equation can be approximated by (in high-energy units): 
E 2 ~ p 2 + m 2 — (E j Eqg)p 2 , from which one can find out that the kinematics of particle- 
production processes (such as the photopion production p + 7 — > p + tt, etc.) will be 
affected - the photopion-production threshold energy gets increased by the deformation 



E > 



(2m p +m 7T )m 7I 
4e 



1 + 



(2m p +m I ) 2 m\ 
64e a E QG 



where E and e are the energies of a 



(m p +m n ) 2 

proton and cosmic microwave background photon, respectively, m p and m n are proton's and 
pion's masses. Similarly, for the electron-positron pair production process 7 + 7 — > e~ + e + 



one obtains the threshold £ > ^ + 
is the electron mass. 



8c*E, 



QG 



where £ is the energy of a traveling photon, m e 



B. Non-perturbative regime 

As was mentioned earlier, if the propagation speed of a particle is very close to c then 
one can not approximate Eq. ( |T5l) by Eq. ( IT6|) . and therefore, the expression ( |T9l) can not 
be valid in general. 

There exist at least three ways of how one can derive here the correct expression for 
velocity as a function of energy. First way is to assume t>2 = v 1 + Av, E% — E\ + AE 
in Eq. ( !T5|) . expand the latter to a linear order w.r.t. Av, replace therein A's by their 
infinitesimal counterparts and integrate in the spirit which led us to Eq. (fT§|) . Second way 

2N 

is to look for an approximate solution of Eq. f lT5|) in a series form Vi/c = ^2 a ra (ej) n , 

n=-2N 

where q = ^Ei/E <^ 1 (E = —£Eqg, as before, £ = ±1), iV > 1 is an approximation 



order, a n are energy-independent constants to be determined. The third method 



271 ] might 



look not as rigorous as the previous two but it is fast and elegant: one should just write Eq. 



f |T5]) in the form J |^| 2 _| = , where its solution becomes obvious. 

All these methods lead to the same result: the desired v (E) is a solution of the algebraic 
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equation a/ (c/v) 2 — 1 = a/x 2 — 1(1 ~~ E/E ), namely 



v/c 



i + (x 2 -i)(i-^ 



1 -1/2 



(21) 



where x — (c/v)e/e =o is the emerging parameter which value can not be determined from 
Eq. ( |T5|) alone. By construction this parameter does not depend on energy of a particle but 
may vary for different kinds of particles. From this expression one can directly computed 
the effective refractive index of the vacuum. In the Cauchy form it can be written as 28] 



n 2 = (c/v^) 2 = 1 + // 7 [1 + M(u)(u /2nc) 2 ] 



(22) 



where /i 7 = x^ ~ 1 an d M(u) = (2-nc/uo) 2 (1 + 2^uo/u) are, respectively, the constant of 
refraction and dispersion coefficient of the physical vacuum, v 7 is the velocity of a photon, u 
is the angular frequency of the electromagnetic wave, u = \E \/h is the natural frequency 
of the vacuum. All this confirms once again that the physical vacuum is the medium with 
non-trivial properties which affects photons and other particles propagating through it, and 
the effects grow along with particles' energies. 

The final dispersion relation can be obtained by subsequent integration, as for Eq. ( 1201) : 

Eq 



P-Po 



2c- s f\i 



T (vW _ E / E o)) ~ arcsinh(y/7) - x^ 



(23) 



where we denoted T(x) = x\/\ + x 2 +arcsinhx, \i = x 2 ~ 1> andpo is the integration constant 
representing the momentum of the background, it is convenient to work in the comoving 
frame of reference where po = 0. It can be convenient also to eliminate x from Eqs. (I2T]) 
and (T23]) to obtain the expression for momentum as a function of energy and velocity: 
E — E 



V 



2c r 



T (17(1 - E/E Q )) - arcsinh(r) - TVT+T 2 



(24) 



where by T = a/ (c/v) 2 — 1 we denoted the inverse Lorentz factor. In fact, this formula is 
the replacement of the relativistic dispersion relation p = Ev/c 2 for ultrarelativistic particles 
with high energies E l-^ol- 

The main feature of the non-perturbative solution is that it indicates the existence of 
the different classes or sectors, depending on the value of x- However, unlike the classical 
relativity, sectors of the "subluminal" (v < c) and "luminal" (v = c) particles are not totally 
disconnected: the propagation speed of the subluminal particles can reach c at finite energy. 
Among other things, this may cause the transluminal phenomena in vacuum discussed below. 
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First mode, called standard or analytic, can be seen when x ^ 0. In this case Eq. (I21"j) 
allows expansion into the Taylor series w.r.t. energy: 

«W y 2 -\E (y 2 - l)(y 2 - -) E 2 

— = l + I ^ +0( E yEl), (25) 

where c x = c/x is the "renormalized" speed of light. This shows that for photons in this 
mode the inverse of x can be interpreted as the "luminal" Mach number and thus x is 
related to the (effective) refractive index of the physical vacuum. To prevent their motion 
from becoming superluminal in this mode, x 2 must be larger than 1 - but not much larger, 
most probably no more than ten percent, as to retain the formal value of c x close to c. This 
makes the dimensionless series coefficients in Eq. (1251) small - in addition to the smallness 
of the ratio E/E Q itself. 

Another mode, called anomalous, is given by the non-analytic branch of the solution (pTj) 
at x — 0: 

,(t) _ c/V2 



l + \^ + O(E 2 /E 2 ) 

4 Eq 



(26) 



\fEjEz 

of course, this expression can be valid when v where energy of a particle can not be 
vanishing. For ultrarelativistic particles there is no obvious boundary condition to rule such 
modes out, therefore, this mode is the (particular example of the) "superluminal" one: it 
is an essentially LIV phenomenon which describes a particle which can propagate in the 
(nontrivial) vacuum with the velocity larger than c. 

Generally speaking, the superluminal modes exist for x 2 < 1, as one can directly see from 
Eq. (1211 . Unlike the tachyons in the classical relativity theory, their energies are real- valued 
and stay finite when v approaches c. If we extend x on to the complex plane (but keeping its 
square real-valued) then the superluminal particles can be further classified depending on 
whether x 2 is a strictly positive number or not. If it is then the minimal allowed energy of 
such particles is zero, otherwise, i.e., when x is imaginary or zero, the energy must be greater 



than some threshold value: E min = E (l — 1/y 1 — x 2 )- Thus, they are not expandable in 
series in the vicinity E = 0. The mode f|26|) is a special case of the second subclass and is a 
kind of the "interface" mode between the subclasses: its E min is zero (similarly to the first 
subclass) but v(E) is not analytic in that point (similarly to the second subclass). 

The common feature of the superluminal particles is for the "sub-Planckian" energies 
(\E/Eq\ < 1) their propagation speed decreases as energy increases - as opposite to the sub- 
luminal mode f )25|) - until it reaches c. This can be explained by when a particle propagates 
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faster than the speed of light in the physical vacuum its interaction with the latter leads 
to the "luminal boom" and appearance of a conical front of the shock waves carrying away 
large amount of energy. In the classical relativity this energy is actually infinite and thus 
the barrier crossing would be forbidden for known particles. This Cherenkov-type radiation 
is mostly electromagnetic but it can lead to creation of other known particles - often with 
very high energies. An interesting question is whether the electromagnetic component of 
this radiation exhibit the anomalous Doppler effect - similar to the one for the superluminal 
(non-point) sources in vacuum which has been predicted even at the level of the classical 
relativity theory 29]. 

From the observational point of view, so far there exist not so much arguments directly 
supporting the existence of luminal booms in the (nontrivial) vacuum in some GRBs and 
AGNs - probably because one can often find more than one way of explaining the majority 
of astrophysical phenomena, especially if a phenomenon is complex and/or the observational 
data lack of a necessary accuracy. Nevertheless, some arguments do exist: for instance, the 
softening of a GRB afterglow is similar to the frequency evolution of a sonic boom: at the 
surface of the shock cone the frequency is very high but rapidly decreases inside, and there is 



even some quantitative similarity 



30( . Also there may appear the phenomenon of mimicking 



the double-lobed radio sources, such as DRAGNs, by such shock waves - by analogy with 
any two sound waves from a supersonic jet initially emitted at different times (and thus 
from different locations) but reached an observer simultaneously thus creating an illusion of 
the doubling of the sound source. In this connection, another interesting question is what 
would be the "luminal" analogues of other trans- and supersonic phenomena, such as the 
Prandtl-Glauert singularity, N- and U-wave, etc. 

To summarize, all this reasoning which led us from Eq. ffT4l) to Eqs. f|T9|) and fT2"Tj) 
indicates that in our theory the dispersion relation for all scales of energy and momentum 
may actually vary depending on a physical situation, and therefore, the complete physical 
picture is still on its way - main reason of which was explained in the paragraph preceding 
Eq. ( |T5l) . In this connection, it would be interesting to find the way of observational testing 
of directly Eq. (TH|) . or Eq. (fT5l) but taking into account cosmological effects, as they are 
the primary predictions of the theory. 
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V. CONCLUSION 



In general, our estimates imply that due to the effects of the nontrivial physical vacuum 
the mean free path of a subluminal high-energy particle, its lifetime in a high-energy 
state and, therefore, travel distance from the source can be significantly larger than one 
would expect from the conventional theory. In fact, using arguments of such kind one can 
show that the deformed dispersion relations above are capable of explaining results of few 
other classes of experiments: observations of cosmic rays above the expected GZK limit, 
studies of the longitudinal development of the air showers produced by ultra-high-energy 
ladronic particles, ATIC observations of the high-energy electrons from an unseen source 
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3lN33|. The trans- and superluminal phenomena are briefly discussed as well. 
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Remarks. After (the initial version of) our paper has been e-printed we found that in 
the very recent article, "A limit on the variation of the speed of light arising from quantum 
gravity effects" by the Fermi LAT and GBM Collaborations, published in Nature 462, pp 331- 
334 (e-printed at larXiv:0908.18 32). the velocity dispersion ffT9|) has been ruled out. For our 
theory it is not a problem though because in that particular case the linear approximation 
described in Sec. IIV Al is hardly valid. For the extremely ultrarelativistic particles, one 
should consider instead the non-perturbatively derived dispersion from Sec. IIVBI eventually 
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leading to Eq. ( )25|) . The latter is not ruled out by Fermi's data - those can only put further 
bounds for the parameter \. 
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